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We search for supersymmetric standard model realizations with extra singlets and extra f/(l) 
using the heterotic string compactification on the Za-u orbifold with two Wilson lines. The effec- 
tive superpotential produced through the vacuum restabilization mechanism is examined for three 
representative Pati-Salam string models obtained in the literature. An automated selection of semi- 
realistic vacua along flat directions in the non-Abelian singlet modes field space is performed by 
requiring the presence of massless pairs of electroweak Higgs bosons having trilinear superpotential 
couplings with massless singlet modes and the decoupling of color triplet exotic modes needed to 
suppress B and L number violating processes. 

P^' PACS numbers: 12.60.-1,12.60. Jv 
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I. INTRODUCTION 



Model building with the heterotic string [1-9] has received a vigorous stimulus thanks to the recent focus [10-26] on 
anisotropic compactifications on the orbifold Zq-jj. Recent reviews are given in [27, 28]. Satisfactory representative 
^j.' , vacuum solutions are also reported for the Z12-1 orbifold [29-33]. The exploration of the vast moduli space of the 
■ Zq-ii orbifold with two authorized Wilson lines [17, 20] was found to give access to a fertile mini-landscape of vacua 
I for minimal supersymmetric standard models descending from 5-d or 6-d grand unified theories. Sampling the still 
. wider space of solutions with three Wilson lines [34] leads to equally hopeful conclusions. 

' One characteristic feature of these top-down constructions is the profusion of J7(l) gauge factors which remain after 
, the first gauge symmetry breaking from the orbifold gauge twist and the presence of Wilson lines. One of them is 
ILj^ ' anomalous, creating a one-loop Fayet-Iliopoulos D-term a bit below the string scale that needs to be cancelled. This 
step usually leaves a lot of freedom due to the large vacuum degeneracy. In the vacuum restabilisation process, the 
U{1) gauge symmetries will get broken and most of the extra modes beyond the MSSM (minimal supersymmetric 
^ ] standard model) matter content, but hopefully not all of them, will decouple with large masses. Except for searches 
of extensions including right-handed neutrino supermultiplets [35, 36], most studies have restricted to solutions in 
which the whole set of Standard Model singlet modes acquire large masses and decouple from the low energy field 
theory. It is natural, however, to ask whether the large moduli space of these compactifications does not include 
regions where extra singlets and extra U{1) symmetries might realize the NMSSM (next-to-minimal supersymmetric 
standard model) or one of the related versions proposed in the literature [37-42] and analyzed over the years in 
phenomenological [43-48] and formal [49-53] studies. (After completion of the present work, there appeared a study 
of NMSSM realizations complying with the severe selection criteria of [17, 20] which reports that solutions exist only 
in the case with three Wilson lines [54].) 

Filling this gap is the purpose of the present work. Rather than pursuing a statistical exploration of the moduli 
space of heterotic string compactifications on the Zq-u orbifold, we restrict consideration in this preliminary study 
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to the three representative 5-d Pati-Salam models, designated in [10] as models Al, A2, B, and apply techniques 
developed in past works. The early studies along these lines focused on flat directions for the non-Abelian singlet 
modes realizing the minimal supersymmetric standard model [55 58] or the Pati-Salam model [59], while subsequent 
studies have explored field directions along scalar modes charged under the non-Abelian gauge group factors [60-63] . 

Having in hand the string massless spectra for a set of representative solutions, we use a Fortran computer program 
to implement an automated search of the flat directions which exhibit a satisfactory supersymmetric effective action. 
Our study is somewhat close in spirit to that of [59]. Concentrating on string theory realizations [64-67] of the 
Pati-Salam model [68] comes with a number of well-known advantages. We shall not detail in the present work the 
Pati-Salam to Standard Model {PS SM) gauge symmetry breaking, but note that the use of Pati-Salam group 
multiplets merely serves a convenient book-keeping purpose, since the representation content with respect to the 
Standard Model gauge group is then uniquely determined. 

The contents of the present work are organized into four sections. Section II discusses general features of the low 
energy theory with Pati-Salam gauge symmetry. A brief review of the supersymmetry preserving flat directions in 
the presence of an anomalous ^7(1)^ gauge symmetry is provided in Appendix A. In order to clarify certain fine 
points that have caused us much confusion, we found useful to include in Appendix B a brief review of the string 
theory construction for the Z^-n orbifold complementing the detailed discussions in [10, 13, 20]. The core part of the 
present work figures in Section III, where we present the results of our searches of singlet extended supersymmetric 
standard models for the Pati-Salam string models A2, B of [10]. A general discussion of results and a summary 
of our conclusions figures in Section IV. 

II. LOW ENERGY THEORY 

A. General features of Pati-Salam string models 

Our interest is on the 4-d orbifold compactifications of the Eg x Ef, hcterotic string satisfying Af = 1 supersymmetry, 
with gauge symmetry group, G x G', including a Pati-Salam observable sector, 

G = G422 X JJ U{l)a, [G422 = SU{4) X SU{2)l x SU{2)n\ (H.l) 

a=l 

and a hidden sector, 

G' = so{mY X su{2y or G" = so{mY, (11.2) 

with the corresponding Abelian factors, U{1)^'', [Na = 5,6]. The gauge group representations consist of three 
chiral generations of bifundamental matter supermultiplets, /j ~ (4,2,1), f^ ~ (4,1,2), a number of electroweak 
Higgs bidoublct supermultiplets, hi ^ (1,2,2), and a number of vector bifundamental supermultiplets, f^, f[ and 
fi'^i fi'^y which accomplish the PS — >■ SM Higgs mechanism for the gauge symmetry breaking, S'[/(4) x SU{2)ii — >• 
SU{3)c X U{1)y, along with a set of non-Abelian singlet modes, (pt ~ (1, 1, 1)(1', 1')- There also occur modes with 
exotic quantum numbers: 

a ~ (6, 1,1), <~ (1,1,2), 4 ~ (1,2,1), gi~(4,l,l), g,~(4,l,l); (II.3) 

hidden sector modes: 

E^~(16M'), E^~(W,1'), 7;'-(10',l'), <~(1',2'); (n.4) 

along with modes charged under the observable and hidden groups: 

/i;'~(l,2,l)(l',2'), h'f ^{1,1,2){1',2'). (II.5) 

Among the set of non-Abelian singlets, (t>i, wc distinguish the modes without and with oscillator excitations, Si and 
Yi. We group the singlets into two subsets, the first including the fields (jii excited along the flat directions with 
large vacuum expectation values (VEVs) of the order of the string or Planck mass scales, rus or M*, and the second 
including the dynamical fields, (j)i and (j)i, with undetermined or vanishing VEVs and masses. 

The PS SM phase transition is assumed to occur at a mass scale nips near the string scale. This hypothesis 
is in harmony with the renormalization group analysis of the gauge coupling constants [10] and with related phe- 
nomenological studies [65-67]. At the PS — > SM gauge symmetry breaking, the above listed multiplets decompose 
into Standard Model group 0321 ~ SU{3) x SU{2)l x U{1)y multiplets as, 

/ ~ (4, 2, 1) ^ g + I = (m d) + (1/ e) - (3, 2) 1 + (1, 2)_ 1 , 
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r - (4, 1, 2)^q' + l' = [d- + (e= v-^') ^ (3, 1)2 + (3, 1)_4 + (1, l)i + (1, l)o, 
h ~ (1, 2, 2) ^ i/„ + iJd ~ (1, 2)i + (1, 2)_ 1 , 

C ~ (6, 1, l)^-g- + g-^ (3, + (3, l)i , 

q^{A,l,l)^g' + E= (3, 1) i + (1, 1)_ i , 

d'~(l,2,l)^(iV' i?')^(l,2)o, 

d'- ~ (1, 1, 2)^E^ + N^-^ (1, + (1, 1) I . (II.6) 

Analogous decompositions hold for the conjugate representations, / ^ (4,2,1), (4,1,2), q ^ (4,1,1). The 

superpotential couplings of Pati-Salam modes of lowest order decompose into couplings between SM modes as 

hkhi = Hu,kHd,i +_Hd,kHu,i, 

QkQi = g'kg'i+EkEi, 

CkCi = gl9t+ glti, 

hfr = Huqu^ + HdqiF + Hdle^, 

Cff = g^qq, 

Cff^ = g'^u'^cf + g''u''e'' + g''dy'=, 

Cff = g^vFd" + g''u''e'' + g''d''v''. (II.7) 



B. Search strategy 

We search vacua of the three chiral generation models Al, A2, B of [10] with the following characteristics at low 
energies: 

• One or several non-Abelian gauge group singlets, (pk, coupled to one or several pairs of electroweak Higgs boson 
bidoublets by trilinear effective superpotential terms, hihj<j)k', 

• Extra U{iy symmetries [48] possibly broken at mass scales much lower than 0(M*); 

• A secluded sector with additional singlets [39, 42] belonging to an approximate flat direction and charged under 
the extra 11(1)' only. The simplest version includes in addition to s three singlets [39, 42] si, S2, S3 belonging to a 
flat direction approximately lifted by a small coupling constant A' in the effective renormalizable superpotential, 

Weff = XHuHdS + A'siS2.s.3. 

The superpotential is decomposed into three main components, Wg, W2, W3, associated to the couplings of the 
non-Abelian singlet flelds, alone or multiplied by two and three fields charged under the Pati-Salam group. Each 
component is built as a power expansion in the of the form 

n n,m i,j 

W2{4'k) = n]{(f'k)fifj + Pi]{4'k)filj + lii;j{(t)k)hihj + Tij{(l)k)CiCj + (Tijqiqj, 

= Cfe[4^(0O/r/; + Ck'i't'^f'if^ + ^i{^i)hf3 + 4{^i)m3 + e^'iMmj] 

+y,^i^i)hkfif^ + Xk'MhkM!, 

(11.8) 

where p\"^\ gj™^ € Z-\., ^ij{p\"^^ +<lj"^^) > n and the summations over fi, (and their complex conjugates) include 
both the matter and Higgs multiplets. The coefficient functions fiij{(f>k), aij{(f>k), Tij{(f)k), ■ ■ ■ are given by infinite 

(n) 

power expansions in the singlet fields of same form as that displayed above for Ws ■ Integrating out the massive 
modes Ci, by means of the classical fields equations, can produce baryon and/or lepton number violating couplings 
represented by local operators of dimension V = A: Iqd'^, u^d'^d'^, lie'', or of dimension D = 5: qqql, u'^d'^d^v^. The 
resulting dangerous local operators, 

Weff = Mt-^i) kf f f f + {lk{r-^)kici) .f.f.r.r + {ck{T'^)kici) ffrr 

-^[l^{(l))qqql + l{(t>)c((l))qqd''u'' + {(t>)u'' d"" d" v"] + ^^c^icpXd'^d'', (II.9) 

may compete with similar couplings already present at the compactification scale. 
Our search strategy consists of four stages. 



4 



• Firstly, we construct (up to some fixed order) the superbasis of holomorphic invariant monomials 

PaW =U^^ = n^r'"^rS [r?, nf, m? e Z+] (11.10) 

i i 

which solve the D flatness conditions, Qa{Pa) = 0, QA{Pa) x sign < 0, where Qa, [a = 1, • • • , A^a — 1] 
denote the anomaly free Abelian charges and Qa the single anomalous Abelian charge with Fayet-Iliopoulos 
parameter ^a defined in Eqs. (A.l) and (B.9). 

• Secondly, we construct the superpotential Wg in the non- Abelian singlet modes (up to some fixed order) satisfying 
the gauge symmetries and the string selection rules in Eqs. (B.IO). 

• Thirdly, we scan the D flat monomials Pa{4>i) and compare each of these in turn with the allowed couplings 
of singlets in Ws{(l>i), in order to determine which of the D flat monomials obey the type A or B F flatness, 
such that there are no allowed couplings in the superpotential Wg with all field factors, or all but a single field 
factor, included in Pa{4'i). (Said otherwise, type A or B lifting of a flat direction Pa{4'i) occurs whenever some 
superpotential monomial with only 1 or no field having zero expectation value is allowed.) Our terminology 
slightly deviates from the original one in [55, 56, 62] which referred to the subset of monomials in Wg allowed 
by the gauge invariance rules alone. 

The conservative approach of exploring the space of supersymmetric vacua consists in restricting to the type B 
directions, on the grounds that these preserve local supersymmetry, thanks to the vanishing scalar potential or 
cosmological constant, < V >= 0, and can be made flat to all orders by imposing a finite number of conditions. 
The existing applications [56, 57, 59, 62] typically stop at the trilinear or fourth orders terms of Ws ('/')• An 
alternative procedure is pursued in [15, 17, 31] which consists in selecting a (typically sizeable) set of singlet 
field directions providing for a satisfactory effective action and solving next the D and F flatness equations out 
to high orders for these singlets. In the present work, we follow an intermediate approach in which we restrict to 
flat directions and superpotential monomials of low orders only but include both type A and B flat directions. 
Since supersymmetry must get broken anyway, we argue that it makes sense to consider the approximate flat 
directions which are lifted by non-renormalizable F terms of reasonably high orders. To motivate our choice 
we note that the condition < Wg >« is automatically satisfied for models preserving an approximate R 
symmetry [70]. We also observe that an initially non F flat direction can be repaired into a flat one by using 
the cancellation with contributions from higher order and/or moduli dependent couplings [71]. One may tie 
the order of the F term lifting operator to the supersymmetry breaking scale F by considering the case of a 
single flat direction lifted by the superpotential, Ws{(j)) = t/'^+'^/Mi'. Assuming that (p picks up from radiative 
corrections a tachyonic soft supersymmetry breaking mass term, V{^) = — m^j^j^, then the scalar potential 
minimization, 

V{ct>) = -m^|<^|2 + i^^|</,P+2|2 ^ ^ v'^\M ^ (m^^'MP)^, (11.11) 

yields the following acceptable range for the breaking scale, |F|^ < (10^ — 10*) GeV, for p = {1 — 3). Larger 
p yield higher scales. It should be noted that the F term lifting by non-renormalizable operators also plays a 
useful role in string theory models realizing gauge symmetry breaking [72, 73] or supersymmetry breaking [74] 
at intermediate scales. 

• Fourthly, we select some supersymmetry preserving flat directions, identifled by sets of cixcited singlet fields 
^i, and evaluate the effective superpotential between the non-singlet modes, while treating the fields VEVs 
as free parameters. Our main focus will be on the bilinear and trilinear couplings of the electroweak Higgs 
modes, hihj{iJ,ij{(f>) + \^j{4))4)k) ■ However, for completeness, we shall also examine the bilinear couplings of the 
modes carrying the bifundamental, sextet and quartet representations, in order to test the decoupling of mirror 
generations and exotic modes and the absence of baryon number violation. 

We pause briefly at this point to clarify some technical issues. The non- Abelian singlet modes VEVs are 0- 
dimensional if they are completely fixed by the gauge constraints, in which case they are roughly set at, (f) = 
0((— ^^/(3a(</>))')- Otherwise, they may be 1-d or higher if they depend on a single or several free complex pa- 
rameters. 

A given holomorphic monomial of the non- Abelian singlets, Pq, = YYi^i 4>i^ , is D-flat if it carries vanishing anomaly 
free charges, Qa{Pa) = Qai't'i)^? = 0, and a finite anomalous charge, Q^(Pq) — QA{4'i)'''f 7^ of opposite 

sign to the charge anomaly, Trace{QA) oc ^a- A superpotential monomial, = IliLi ^i' ' allowed if it carries 
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vanishing anomaly free and anomalous charges, J2^=iQa{^i)^i — 0' J27=iQA{(f>i)sf = 0, and satisfies the string 
selection rules listed in Eq. (B.IO). 

We define the order n of a holomorphic D flat or superpotential monomial, Pa{4>i) or Wx{(f>k,<Pi), by the number 
of field factors this contains, independently of the values of the fields positive integer exponents, rf, s^. The number 
of inequivalent monomial solutions rises fast with the order n considered, reaching orders of hundreds already at the 
relatively low order, n = 4. The solutions will be ordered according to increasing values of the number of field factors, 
n = 1, 2, • • •, while allowing the power exponents to range over the discrete set of values, rf , sf e [0, 1, 2]. Note 
that with this definition, the order n of a monomial is always bounded by its effective order, Y^^=i or *f ■ To 

avoid dealing with an ovcirwhcilming number of solutions, we generally restrict to orders n < 4. 

The unbroken (anomaly free) Abelian symmetries U{l)x along a flat direction are identifled by the kernel of the 
charge matrix, Aai = Qa{4'i), [a = I, iV^ — 1] where ranges over the modes excited along the flat direction and 
the index label a ranges over the anomaly free Abelian symmetries. In practice, we proceed by evaluating the left 
eigenvectors with zero eigenvalues of the matrix, Aai, [a, « = 1, • ' ' j -^a — 1] namely, x^A = 0, for some selected subset 
of the and next retain among the Abelian charges, = J2a ^aQa, those yielding zero charges, Qx{(t>i) = 0, for 
the full set of (jii. 

The effective trilinear interactions descending from non-renormalizable operators of order n + 3 are assigned string 
tree level coupling constants given by the approximate formula [57] 



' (ttMO" 



ijk 



(11.12) 



where C„ arc constant coefficients of 0(1) and /„ = jJXi^i(Pzif{zi,z.i) are integrals over the location of vertex 
operators on the world sheet sphere surface. The integrals for the four and five point string amplitudes evaluate 
numerically to [75-77]: /i ~ 70, I2 — 400. The string coupling constant is commonly set at, Qs = Qx/V^ — 1/2. For 
the case of a 0-d flat direction, assigning the singlet fields the common VEV <j) determined by the Fayet-Iliopoulos 
term, one infers the following formula for the coefficients 



n/2 

|,| , U y 9sM, I TracejQA) \ ' g^+'ajn ( TracejQA) ] 

l"^!-* QA{ct>)J ~ ^ \ 192QA{4>)VkIJ ''^'^ - (7^2^A92)" 1, Q^(<^)^/fcI^' ' ^ ' ' 

With increasing order of the non-renormalizable operators, we may thus anticipate a strong or mild suppression 
depending on the unknown extrapolation of /„ for large n and the model dependent size of the anomalous charges. 
Using, say, /„ = 10'', C„ — 1, and dropping the factor depending on Qa by setting, ^^^^^^^2^ gives, 

eg ~ 10~^" <?l^"C„/„ « j^Q-2n+3^2"+i. A similar formula holds for the /i-term coefficient of the electroweak Higgs 
bosons bilinear coupling, c^^-hijij, by substituting, c-"^ = /i^"^-^*- Bounding the coupling constant by the requisite 
ratio between Planck and effective supersymmetry breaking mass scales, At("V-^* < TeV/M* ~ IQ-^^^ yields, n > 8, 
which shows that non-renormalizable effective operators starting from 0{(j)''^^)hihj can be tolerated. Considering 
instead the tentative extrapolation, /„ ~ 10'?, [C„ = 1], and the numerical estimates for the anomalous charge 
appropriate to model Al (to be detailed below) 

1 

\1\ 1 / ^„„„„^^ .n\ 2 

0.23, (11.14) 



then the corresponding condition on the /i-term coupling, p!^"^^ jM^ — 5 10 " < 10 , gives the stronger bound on 
the order of non-renormalizable operators, n > 14. Note that lowering 0j should weaken the bounds on n. 




III. SINGLET EXTENDED SUPERSYMMETRIC STANDARD MODELS 

The Pati-Salam models Al, A2, B of [10] for the orbifold Zg-// of with two Wilson lines only (since W'^ = 0), 
are described by the 16-d shift vectors for the gauge twist and Wilson lines of the Eg x Eg group lattice: 

Model Al: y=^(2,2,2,05)(l,l,0«), = ^(3, 0^ 33)(0«), = ^(2, -2, 0<^)(0^ 4, 0^). 

D DO 
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Model A2: F = ^(2^, O^)^^, 0«), = ^(3, 0^ 33)(0«), W^s = ^(4,2,-2,0^)(0,4,22,0^). 

ODD 

Models: y= ^(4,l2,0^)(22,06), W2 = ^(0^ 3^ 03)(3, 0, 3, 0^), W3 = ^(0,2,-2, 0^)(0^ 4, 0^). (III.l) 

DO D 

The 4-d gauge groups for the above three models are: 

Model Al : S'C/(4) x SU{2)i x 5J7(2)2 x 1/(1)^ x (S'O(IO)' x 5'C/(2)'). (III.2) 
Model A2 : 5'C/(4) x SU{2)i x 5?7(2)2 x U{lf x (S'O(IO)' x SU{2Y). (III.3) 
Model B : S'C/(4) x S't/(2)i x SU{2)2 x x S'O(IO)'. (III.4) 

Our results for the massless string spectra of models Al, A2 and B match to those of [10], except for certain 
assignments of oscillator numbers, Njj, hence of Ri charges. The numerical data for the linear and cubic traces 

of the anomalous U{1)a gauge charge and the Fayct-Iliopoulos term, ^a/ (SsMi,)'^ = Trace(Q/i)/(1927r^-y/A!^) are 
displayed in the following table. This also includes results for the total number of massless string states and the 
degeneracies of the various modes. 



Model 


Tr{QA) 


Tr(Q\) 


TriQAQl)/2 


kA U/{gsM,f 






hi 


C-i fi 


.h 




/f 


{(U-(li) 


Al 


-2 


-0.00437 


-0.16667 


0.01750 -0.8 X 10-2 


78 


34 


7 


9 3 





5 


2 


4 


A2 


-131.4 


-393.90 


-10.94 


24 -1.4 X 10-2 


77 


36 


3 


3 3 





4 


1 


6 


B 


13.66 


0.312 


1.139 


0.182 1.7 X 10-2 


84 


48 


6 


3 6 


3 


6 


3 


4 










A. Model A\ 
















The 


massless 


string 


spectrum 


of model Al is 


displayed 


in 


Table 


I. 




The 


set 



holomorphic invariant monomials of order n < 2 consist of eleven solutions. Pa = 
[FiiSg, 53^6, YiaY^s. YiAY^^i, ^isifg, Yi^Y^^, Y^^, FigFig, Y^^Y^^. Y^^Y^^, Y^^]. None of these D flat di- 
rections is found to be lifted by the superpotential couplings in Wg at orders n < 4. Rather than analyzing each 
monomial individually, we combine these multiplicatively into what we term as the 'composite' D-flat direction, 
•p(2) — . . .Pjj^^ which then includes the eight modes, Yn, 5*8, 6*3, 5*6, Yi4^, Yis, I19, Yib- The set of excluded 
couplings along the composite direction is clearly the union of all the sets excluded by the individual couplings of 
order n. All these directions are 0-d ones, hence with a common VEV given by, |(?!)|/(g(sM*) ~ 0.23. The effective 
superpotential components obtained by assigning all the flelds along V^"^^ a common value (f>i = (f> read at order n < 4 

10'S'l5- 

Wh = hihi[^''Y„ + cf^^S^] + h2h2[(f>^Sio] + h^hsiSi] + hih2[(f>^ Sw] + hih^cf^^ + Sg] + /l2/l3[0' + </'^5'9]. 

Wc = CiCi[0-Vio] + C2C2[<^^5io] + CaCsiSi] + CiC4<f>''Yr] + CiC2[0] + dCai^'^Ya] + Cidlct^YaYQ] 
+C2C3[^' + (jy^So] + C2C4q)^Yg] + C3C4Y5S11]. 

Wg = qiqi[0] + qiq2[(l)^] + q2qi[0] + q2q2[Si4]- 

Wf = mm + /f/2[52] + f^fm + f^m^'Ys] + f0nsi2] + mm- 

Whffc = /ii/i/f [1] + hrhmm + /ii/i/3i</''F9] + /11/2/f [0] + /11/2/2I0] + /11/2/3I0] 

+/i2/i/ii'^'] + h2hm''Yg] + h2hm<i^^Yg] + h2.f2.fm + h2f2.fm] + ^^2/2/310] 

+/i3/i/f [0] + /13/1/2IO] + h^hfllY^] + h3hr2[Yi] + /i3/2,/f [Fi +Y2 + + h3,f2m^\ (III.5) 

The results for Wh through Whff are illustrative ones in the sense that they include only a subset of the couplings 
at most linear or quadratic in the singlets which are not part of the flat direction. For convenience, we have listed 
the results for the bilinear couplings along P^^) ^\iq second column of Table II. The flavor bases for the massless 
fermions consist of the 3 modes /i,2 and the 3 linear combinations of the 5 modes /f 2 3 which are left unpaired. The 
fermion mass matrices are sums of 3 x 5 matrices < hk > /i,2/i 2 3) [f^ = 1) 2, 3] involving the combinations of 7 light 
bidoublets which acquire flnite VEVs at the clectroweak transition. The presence of an unsuppresscd coupling hi /i/f 
suggests that hi should be the dominant component of the light electroweak doublets and /i, should be associated 
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TABLE I: Spectrum of string massless modes for the Pati-Salam model Al of [10] with gauge group, G422 x (SO(IO)' x 5?7(2)') x 

U{lf. The pair of group factors 5?7(2)i,2 identify with SU{2)l,r- The non-Abelian singlet modes ~ (1, 1, 1)(1, 1)' without 
and with oscillator excitations arc denoted by Si and Yi. The non-singlet modes include the doublets, d\ ~ (1, 2, 1)(1, 1)', dl ~ 
(1,1,2)(1,1)', d'i ~ (1,1,1)(1',2'); the bidoublcts, ~ (1,2,2_)(1, 1)'; the quartets, q, ~ (4, 1, 1)', qi ~ (4, 1, 1)(1^1)'; 
the sextets, d ~ (6, 1, 1)(1, 1)'; the bifundamentals, /f ~ (4, 1, 2)(1, 1)', /f ~ (4, 1, 2)(1, 1)', _/i - (4, 2, 1)CL^1)', fi ~ 
(4,2, 1)(1, 1)', ; the decuplets, T/ ~ (1, 1, 1)(10', 1'); and the scxtuplcts, - (1, 1, 1)(16', 1'), S- ~ (1, 1, 1)(16', 1'). The 
first column lists the twisted sectors label g; the second column lists the modes name; the third column consists of three sub- 
columns which list the twisted subsector specified by the excited Wilson lines, with A'32 = 1, [g = 0]; N32 = 3n2 +713 + 1, [g = 
1]; A^32 = na + 1, [g — 2,4]; and Ns2 = n2 + 1, [g ~ 3], the complex phase 'y{g) describing the orbifold twist action 
on the fixed points in the Ti plane and the modes multiplicity (degeneracy) V; the third column consists of three sub- 
columns which list the shifted -momenta specified by the charges, 6R' , [I = 1,2,3]; and the fourth column consists of 
five subcolumns which list the Abelian gauge charges in the rotated basis starting with the anomaly free charges Qa and 
ending with the single anomalous charge Qa- In order to quote integer charges, the Abelian charges have been rescaled as, 
(5i/0.1667, Q2/O.O2O83 , (53/0.0104, 4 X Q4/O.OO6944, (5a/0.008333. (Note that the first column entry for g is omitted in the 
right-hand part of the table.) 
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to the third generation fermions. The degenerate pairs of modes /2 and /| must then be associated to the first 

two generations. As long as the Dihedral D4 discrete symmetry of [10, 86] is unbroken, the fermions mass matrices 

/2/f would include pairs of identical rows, which means that their rank is < 2. This seemingly good starting point 

motivates one to search for mechanisms breaking the D4 symmetry. As discussed in [10, 86], the non-renormalizable 

operators, hkfifjW^^i^^i SkYiO'r^^,, involving the composite singlet fields, Of^ = ffjj, produce contributions at 

the PS —7- SM transition which can improve predictions for the fermions mass matrices. We note that an alternative 

way to break the D4 symmetry is by switching on the third Wilson line, n2W2- 

(2) . . 

That no constant or linear terms appears in Ws indicates that all the n = 2 flat directions are indeed unlifted. 
The decoupled massive singlets are, 5*1^10,11,13,15 and is, 9. The //-term bilinear coupling in Wh is a rank 2 matrix with 
the massless mode, (/i2 — (p^hi). However, the suppressed component is along hi, which is the preferred mode whose 
VEV supposedly dominates the fermion masses. The trilinear couplings W = ((/)^/i2^2 + 0^/'-i^2)'S'io are seen to select 
the massive singlet ^lo, which figures among the modes which acquire large masses along V^^K However, it should 
be noted that 5*10 remains massless for a subset of the n = 2 directions Pa- The above result for Wf shows that no 
mass pairings among the matter modes are allowed. The mass matrices Mfjfifj have the single non- vanishing entry, 

M/j =< /ii >, if we omit the massive /13 mode. Examination of the exotic modes mass couplings, shows that qi, q2 

and C2, C3 pair up, while the pairs (72, 9i and Ci, C4 remain massless. We also find that except for the coupling, 
Ci/2/2 ['^■^]> all other trilinear couplings, Ci/f/J, Cifffj in Wcff are suppressed. 

The 4x8 matrix of the Abelian charges for modes along V^^^ is found to have column rank 2. The existence 
of a secluded sector requires finding three singlets, 0/, cfij, 4>k, charged under the unbroken U{V)x with a mildly 
suppressed trilinear coupling, c{4>k)4'i4'j4'K- A wide range of choices is clearly possible given the sizeable number of 
singlets orthogonal to the moduli space of V^'^\ and the abundant number of monomial solutions in Wg. An inspection 
of the allowed couplings lets us select the following effective cubic or quartic order couplings yil2>5'i2, Y^SrS 12 < 
Si4 >, < Y4 > Yi2yi6'S'i2, YiQSrSi2 < Si3 >, that can give rise to a secluded sector with three massless singlets 
coupled by mildly suppressed trilinear interactions. A systematic search does not appear warranted at this stage. 

The special role of Sio motivates us to examine its low order self-couplings, which are given for n < 4 by VFs(5'io) = 
Sio[0] + SiQ[(j)^S7] + SiqIcP'^S^]. These resuls indicate that a cubic self- coupling, S'lg^'* < >, cannot arise without 
an unacceptably large mass term, Sf^cj)^ < Sr >. Since this correlation between the quadratic and cubic couplings of 
^lo also holds for the individual flat directions, we conclude that the occurrence of a trilinear coupling, hihjSio, with 
a finite cubic self-coupling is not the favored option. 

We also found useful to study a class of solutions on an individual basis. We present in Table II a subset of the 
bilinear couplings for four randomly selected F flat monomial solutions of order n = 4. The monomials Pj^\ Pni^ Pjv 
are seen to allow a pair of massless bidoublets, hi, /12, with preferred trilinear coupling, Wh = hih2SiQ. Along Pi'f, 

the diagonal trilinear coupling, Wh = <?i^/i2/i2'S'i3, is also present. The direction Pj^^ has all three bidoublets massless. 
However, some subsets of the exotic modes always remain massless. 

To complete the present study, we have scanned the 70 and 48 flat monomial solutions of orders n = 3 and n = 4 in 
search of the bidoublets couplings of form, Wh = hihj{iJ,ij + A^^-^*'). The global results combining the contributions 
from the various individual flat directions read: 

.n = 3: Wh = hihi[0] + h2h2[^^Sio] + h3h3[^Sii]+hih2[^''Si3]+hih3[^''Yi5] + h2h3[(p^ + (l>''Yi5]. 
.n = 4: Wh = hihi[0] + h2h2[(l>''Si2,i3] + h3h3[(t> + (l)^Yi2,i4] + hih2[Sw + cl>^Ye^ii] 

+hih3[S, + ^'Ye] + /i2/i3 [0^^2,14]. (ni.6) 

Happily, it appears that the cases with no trilinear couplings, ^ 0, Xfj = 0, are outnumbered by those with no 
bilinear couplings, jj-ij = 0, A^^ ^ 0. The results for the directions n = 3 favor a scenario with a single light bidoublet 
hi having trilinear couplings, (jP'hih2Si3 or (jPhihsYi^. The results for the directions n = 4 favor three light bidoublets 
with canonical and non-canonical type trilinear couplings. 

B. Model A2 

The mass spectrum of model A2 is displayed in Table HI. The D flatness conditions are solved 
at order n = 2 by the single monomial, YgSg, and at order n = 3 by the 7 monomials, 
SiSiSi5, YgSiSn, SsYieYis, SsYirYis, 53^17^9, 53^6^9, SgS'ioS'is- The number of solutions at order n = 4 
exceeds O(IOO). The combined composite flat direction V^"^^ x V^^^ for n < 3 includes the eleven modes, 
Ig, Sg, Si, Si5, Sii, S3, Y16, Yis, Yig, Yn, Siq. This remains unlifted by all the singlet couplings in Wg of 
order n < 4. However, no unbroken U{1) charges survive if the above eleven modes simultaneously acquire finite 
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TABLE II: Superpotential couplings for model Al of bilinear order in hi, Ci,qi, qi and orders n < 4 in the singlets. The 

column entries refer to the composite flat direction "P'^^ and the four randomly selected individual flat directions: Pj^^^ — 

9 (A) 9 (4^ 9 9 (4:) 9 

YiiYieYii), P}j' = SiS3Yi6Y{g, Pf/j = SiS^S4Si, P^y = 5'i54yi8. Empty entries correspond to cases where no coupling is 
present up to order n = 4. 
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VEVs. The effective superpotential component Wg and those of bilinear and trilinear couplings obtained by assigning 
a common VEV (p to the modes c^i € V'^^^ x T''^^ are given by 

Ws = [(t>^S<,] + [<t>SiS<, + cp\Y^ + Yi+Yi) + cp''YeSir + (l)''Si3Sn + <t>''yioSi3]. 

Wh = /li/li[<^=^ + 0^y5] + /i2/l2[O] + /ll/l2[<^'i3]- 

Wc = CiCi[<p''Y3] + C2C2[<l>Ye] + CiC2['p''Y3]. 

Wg = qiqi[4>^Y^ + </.3y2] + q2q2WY^] + q^q^[4>ylY^ + (/.F2I7] 

+qiq2[<l) + ^Sq + ^^Fis] + qiq3[(t>'^Yi] + 92^3 [0] + 9192 [0] + qiq^ifS^^Sn] + «2Q3[</'^yiy7]- 

Wf = /f/r[</'r3(ri2 + ri3 + 5i3)] + /27f[o]. 

T^W= = /ii/i/f [^3^4] + /ii/2/210] + /ii/i/|[0] + /ii/2/f [0] + /12/1/f [0] + h2f2mYw] + h2hfm 

+ /i2/2/f[0]. 

Wcff = c^fum + Ci/27|[0'>i4] + 02^1^12] + ^2/2^2 [<^'n]. (in.7) 

We sec that the singlet superpotential Wg has no constant term but includes a harmless tadpole for the massive field, 
^g. The massive singlets are 5'i,9,i3,i7 and 1^1,2,4,6,10- The Higgs bidoublet couplings in Wh show that h\ decouples, 
while /i2 is left massless but lacks trilinear couplings to singlets. Indeed, the quartic coupling, (^yio5i3 /12/12, cannot 
generate an effective trilinear coupling at lower scales, since both of the singlets Yio, 5*13 are massive. No large mass 
pairing takes place between the conjugate bifundamental modes. The suppressed couplings in Whj preclude the 
possibility of discriminating between the light and heavy flavors. Some zero entries may be possibly lifted by assuming 
a finite VEV for the singlet Yio. The sextet and quartet modes, except the pair, 91^2, remain massless. The couplings 
in Wcff give no contributions to baryon number violating operators. 

We have also scanned the individual D flat directions of orders n = 2, 3 and n = 4 in search of the bidoublet modes 
couplings. The full set of 15 monomials at orders n = 2, 3 are F flat. Of the 100 monomials at order n = 4, about 
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TABLE III: String spectrum of massless modes for model A2 of [10] of gauge group, G422 x (50(10)' x 5[/(2)') x U{lf, 
with the pair of group factors, S(7(2)i,2, identified with SU{2)l,r- Same notational conventions as in Table I are 
used. The Abehan charges (given by integers except for Q4) are defined in the rotated charge basis by the rescaUngs, 
Qi/0.1667, Q2/O.I667 , 4 Q3/O.OO26O, Q4/0.0007667, 18 Qa/0.1244. 
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80 are lifted by F-terms. The results obtained by combining the contributions from the various flat directions are 



We see again a clear trend in favor of several light bidoublets coupled by trilinear terms to massless singlets. 



The massless spectrum for model B is displayed in Table IV. Wc find no holomorphic invariant monomi- 
als at order n = 2. At order n = 3, there appears 16 solutions of which a few representative monomials are: 
Pa = [S11S1SY25, S'ii5'i8l26, 5'iiy23'S'20, 5'ii i^24'S'2o] . The composite direction, 'P'^^^ consists of the eiglit modes, 
^ii. Sis, Y2Q, Y23, S20- Y24. 5i5. No Abehan charges remain unbroken along "P^^^ Assigning a common VEV 
(j) to these modes yields the reduced superpotential 



Ws = cj^YsSs + CPY4S5 + (pSsYr + cj,S5Y5 + cj^SeYio + cj>S8Ys. 

Wh = hihi[Sie] + h2h2[0] + /is/isM + /11/12M + /ii/i3[52] + /12/13M. 

Wc = C'iCi[0] + C2C2[Sie] + CiC2[0]. 

Wq = qiQi [(j)] + q2 52 [0] + 5i 92 [4>\ + 92 Si [0] . 

Wf = /l/l[0] + /i/2[0] + /2/l[0] + f2f2m + cj^Ys] + /i7"f [^l(yi9 + Y20)] + f^fm 

+f2msi] + fUi[SiY2i] + m[SiYi9] + fU2[4'^yiY7]. 

Whffc = hihmYiQ]+hihf^[Yi5 + Yi,Yi^] + hi,hf^[^Yi2 + <l>Yu] 

+hif2fm + hif2fm + /ii/2/i[o] + /12/1/f [0] + h2hfm + h2fifm 

+h2f2fm + /12/2/2IO] + /12/2/3IO] + /i3/i/f ['^'^2 + 4>Yu] + h3hf^[Yi5] + hahfilYu] 

+hsf2fm + hsf2fm + hsf2fm- 

Wcff = Ci/i/i[0] + Ci/2/2[0] + CiA7f [Fig + 1-20] + Cif^mYis] + Cif^mcpYi2 + <pYi4 

+C2/i/i[o] + C2/2/2[o] + C2fum + C2fum + C2fifm- (ni.9) 



The identical couplings for hi and C2 is a consequence of the fact that these modes have same quantum numbers. The 
singlets superpotential Wg contains no constant or linear terms, while bilinear terms appear for the singlets £'3,5,6,8 
and 13,4,5,7,8,10, which thus pick up large masses. The bidoublets couplings in Wh show that the three sets of hi modes 
remain massless, and have trilinear couplings {hhS2 + hhSio) involving the massless singlets 5*2 and Siq. No mass 
pairings arise between the conjugate bifundamental modes that would remove some of the fi, ff modes or decouple 
some of the mirror fermions. fi and /f, ff. The bidoublet-fermions couplings in Whff" give vanishing mass 
matrices for the matter fermions. The flavor structure is problematic since it gives no hint on how to discriminate 
between the heavy and light flavors or between the matter and Higgs modes, independently of the dominant linear 
combinations of the bidoublets. The results for Wc and Wq show that all the sextet modes, C,, and the single pair 
of quartet modes, q2, qi, remain unpaired. 

The special role of the singlet Siq motivates us in studying its self-couplings at orders n < 4. The results for 
7^(3) . \Yg(^SiQ) ~ 5'i6[0] + S'ff5[0] + S'j'gp], indicate that the mode 5*16 is light but has a suppressed cubic coupling. 
The results for V^'^\ W,{Sie) = Sie[0] + Sf^lcp^Yis + <P^Yii] + Sf^icp'^Ys + 04 Yg] would also allow for a light mode Sw 
provided that the massive singlets 13,8,11 acquire small VEVs. 

Results for seven randomly selected individual flat directions are displayed in Table V. For all monomials other 
than Pjv, we see that all the hi remain light and have trilinear couplings to singlets, but that the sextets and a subset 
of the quartets fail to decouple. Wc have also performed global type scans of the D flat monomial solutions restricted 
to the bidoublets couplings. Of the 16 and 100 solutions at orders n = 3 and 4, we find and 10 monomials lifted by 
F-terms. The bidoublets couplings combining the contributions from the various flat directions are given by 



• n = 2, 3 : Wh ^ hihi[Yis] + /i2/i2[0] + hih2[0]. 

. n = 4 : Wh=^ hihi[4>''Yis,i9] + /i2/i2[<^^ll4,i5 + ^^£13] + hih2[4>^Si + ^"^Y^ + .^^^^J. 



(III.8) 



C. 



Model B 



• n = 3: Wh = hihi[Sw]+h2h2[Q\+hzhi[Q\ + hih2[Q\ + hih^[S2] + h2hz[Q\. 

• n = 4: Wh = hihi[Sie + 4>Yi9 + (l>Y2o] + h2h2[0] + h3h3[<j)Yi5 + (j)S2] + hih2[0] 

+hih3[(f> + 4>S2 + (f>Sie] + h2h3[(f>^Yi3 + cj>''S^ + ^^Fi]. 



(III. 10) 
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TABLE IV: Massless string modes for model B of [10] with gauge group, G422 x {30(10)') x U{1)^. We 
use same notational conventions as in Tabic I, except that the pair of group factor 5'[/(2)i,2 identify here 
with SU{2)r^l- The quoted integer rescaled charges are related to the initial ones by the rescalings, 
4 Qi/0.005208, 25(52/0.002042, 12 Q3/O.OO96I5, 9 Q4/O.OO347, (55/0.000672, gA/0.017796. 
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TABLE V: Bilinear superpotential couplings for model B in the modes hi, Ci,qi, qi of order, n < 4 in the singlets. The 
column entries refer to the three fiat directions of order n = 3 : Pj = S11S18V25, Pii = 511123520, Piii = S'ir,l23l2(i- and the 
four randomly selected flat directions of order n = A : Piv = S2Y6SfiS20, Pv = YiYi^Si7Y26, Pvi = YtiSi^Yi^jYzQ, Pvii = 
IgySisyisVis- Empty entries correspond to cases where no coupling is present up to order n = 4. 



w 


Pi 


Pii 


Pin 




Piv 


Pv 


Pvi 


Pvii 




S16 


SlQ 


S16 




S16 + <l>Yl9,20 


S16 




SlQ 


h2h2 




























0Vl5 


<i>S2 + 0^Sl6 


<j>S2 + <l>^Sie 




hi/i2 


















hihs 


S2 


S2 


S2 




<P 


S2 + 4>Sl6 


S2 + 0Sl6 


S2 


'12/13 










<l>^Yi3 


<t>'^Yi + 0=56 






CiCi 


















C2 C2 


S16 


SlQ 






S16 + <t'Yl9,20 




5^16 


SlG 


C1C2 


















qiqi 











Y25 


























9192 


Y23 








4> + </>^yi2 


^25+0^^11 +0^5io 125 


+ 0^Fii +<^^Si9 Y25 + 


?i^F22 +0*Y-19,2O 


92 91 










Y23 


Y23 


Y23 


<^ + <A^y2i 



These results again favor the scenario in which light bidoublets interact by trilinear couplings to singlets. A scan of 
the flat directions reveals the existence of a large majority of individual monomial directions along which the mode 
^le is light but with a suppressed cubic self-coupling. 

IV. DISCUSSION AND CONCLUSIONS 

Let us first state some general features of our results. The top-down constniction for the Zq-jj orbifold has a large 
number of flat directions and a rich structure of superpotential couplings. This contrasts with the situation prevailing 
for the Z3 orbifold models [1, 7, 35], but is in harmony with that in non-prime orbifolds or free fermions. Comparing 
with the string spectra obtained in intersecting brane models [69, 78] would not be very teaching because the existing 
type I string models are more akin to bottom-up type constructions. 

We have carried the F flatness test indiscriminately for both types A and B flat directions, restricting to the 
superpotential monomials with at most four distinct singlet field factors, Ws{(f>i) = YYi=i 4'T ■< — ^' — •^]- These 
contain a subset of the monomials of absolute order < 2n = 8. It is intuitively clear that the F flatness condition 
is more severe for the lower order superpotential monomials, since these are more likely to contain all (or all but 
one) of the fields excited along the fiat direction. Upon pushing the F fiatness test in model B to the order n = 5 
of Wg, we find that out of 100 order n = 4 flat monomials, 34 are lifted by the 3 monomials of order n = 3 in 
Ws = SiYiQ, 1455125, S5Y5Y2S. Including all the monomials in Wg up to order 8, selects, 0, 6, 40 flat monomials of 
orders n = 2, 3, 4. 

To determine which of the singlets decouple and which remain massless, we have derived the mass matrices for the 
Si and Yi by scanning over each individual flat direction \ip to absolute order 8. We found that one cannot decouple 
all the singlets at the string scale by these means in any of the models. However, taking the combinations of all the 
singlets of a given order, we found that all singlets of models Al and A2 can be decoupled if the order is large enough, 
n > 5 and 7 respectively, and the number of fields involved is itself large enough. One might prefer vacua defined by 
smaller sets of fields, and indeed our aim here is not to decouple all the singlets but leave a few of them massless in 
order to obtain an NMSSM like model. We thus prove that there is room for decoupling all but a few relevant modes, 
should one perform a thorough scan over the fiat directions. As far as model B is concerned, the number of singlet 
fields is too high to decouple all of them at the order considered, and making this model viable on a phenomenological 
basis would require giving a lot of singlet modes a mass at the level of the Pati-Salam breaking down to the Standard 
Model gauge group. 

Our study of the fermion mass generation has been rather sketchy. The preferred candidate for the electroweak 
Higgs bidoublet is obvious only in model Al. To discriminate between the heavy and light fiavors in models A2 
and B requires a closer analysis of the couplings of bidoublets and bifundamentals. Such a task would be warranted 
once one has really in hand a benchmark type model. All three models satisfy at the string mass scale a D4 family 
symmetry with a number of bifundamental modes / or f from the twisted sectors Ti or T3 transforming as doublets. 
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An acceptable description of the fermions flavor structure can be achieved only if the D4 symmetry is broken. On side 
of the promising mechanism using the condensation at the Pati-Salam breaking scale of the composite singlet fields, 
O'lj = f[fj, stringy mechanisms can be envisionned to lift the degeneracy of the modes in the Ti^3 sectors. One could 
use blow up submanifolds of different sizes at the two fixed points, or consider the so far lightly explored models with 
three Wilson lines [34]. 

The existence of a light pair of Higgs bosons seems to be correlated with the presence of massless exotic color 
triplets descending from sextet and bifundamental modes, Ci, f^'^. This is unavoidable in models Al and B because 
of certain hi and Ci having identical quantum numbers. A single pair of g^, qi modes always fails to decouple. Since 
the exotic modes are vector like, they could pair off or couple to lighter modes into which they decay, through lower 
scale physics, including the PS — >■ SM transition. We have not examined here the decoupling of the weak doublet and 
singlet exotic modes, S = {N\E''), dT = {E^,N^), since these represent perhaps a lesser threat on the low energy 
theory. 

Our main purpose in this study was to establish an existence proof, based on three representative string models, 
of supersymmetric models where suppressed bilinear couplings of the elcctroweak Higgs doublets coexist with unsup- 
pressed trilinear couplings to singlets, jjLhkhi + Xhkhifj), [ji < 0{(j)^),]. Rather than scanning over the flat directions, we 
have followed an approximate procedure making use of what we term as composite flat directions. This is admittedly 
open to criticism, especially at high orders, because of the large number of excited singlets. However, the qualitative 
orientation this approach gives does not seem invalidated by selecting randomly, or scanning in a global way, a subset 
of individual flat directions. We frequently encounter canonical and non-canonical type trilinear couplings, Xfjhihj4>k, 
however, without any obvious correlations between these bidoublets and those dominating the fermion-Higgs Yukawa 
couplings. 

Additional 11(1)' gauge symmetries are generally present in the individual flat directions of order n < 4, but they 

are absent in most composite directions. It is difficult to decide if a gaugino supersymmctry breaking mediation or 
a secluded sector scenario is favored since our searches give no clue on the size of the Z' boson mass scale and the 
Z — Z' mixing angle. 

In summary, the search of models with extra singlets is considerably facilitated in the Zq-jj orbifold by the rich 
mini-landscape of vacua. We find a clear preference towards an active role for extra singlets, but the presence of 
sizeable trilinear couplings hhcj) along with strongly suppressed bilinear couplings hh is not systematic. The study of 
solutions was made tractable thanks to the restriction to low orders, n < 4. Beyond this order one must consider 
representative samples of the flat directions. However, the composite flat directions of low orders, n < 4, seem to 
capture general features not invalidated by global scans of the individual directions, as confirmed by the similar 
conclusions for models Al and B. The strong 0{(j)^) suppression required for the /z couplings suggests, however, that 
definitive conclusions could only be made until the searches include higher order couplings. Weak and strong points 
are present in all three models, with no model faring best on all issues. 
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Appendix A: Flat directions in heterotic string compact ificat ion 

Given a gauge theory with the Lie group Ga and superpotential W, which includes the set of massless chiral 
supermultiplets, 4>i, [i = 1,---,A''] then the D flatness conditions are given by: = = {4>i{T'^)ij4'j)j where 
are the Lie algebra generators, and the F flatness conditions are given by: Wi = = 0. In the scalar fields 
(j)i vector space, the D fiat directions are parameterized by the gauge invariant monomials, Pa{(j>i), satisfying the 
equations, Pa^i = §^ = Cac/)*, for suitably chosen complex constants, Ca, with the (f)i substituted by their (constant) 
VEV [1, 80]. (This statement is easily checked by writing the Ga gauge invariancc condition as, = SaP,, = 
J2i T^^afpi = Ca J2i <f>i{Ta4>)i-) For an anomalous gauge symmetry [81-83] U{1)a with the charges, Q\ = QA(<Pi}, 
and a finite charge anomaly implying the presence of a Fayet-Iliopoulos term with coupling constant S,a, the D-flatness 
condition is modified to 

= DA = Y,WA^i+U, 



«2 A/f 2 1 

U = e'^M'j^s = f^^Tr(Q^), kA = 2^(Qi)^ 



(A.1) 
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where gs =< e*^ > is the string couphng constant, with (j> the 10-d dilaton field, and Q^, [/ = 1, • • • , 16] are the 
orthogonal basis components of the anomalous charge generator in the Eg x group weight lattice. For the set of 
gauge group factors, G = (Ha ^a) x U{1)a, the D flatness conditions is solved by considering holomorphic invariant 
monomials, Pa{(f>i), uncharged with respect to the anomaly free gauge factors Ga, but with finite anomalous U{1)a 
charge, QAiPa) 7^ 0, of opposite sign to oc Trace{QA)- In the supergravity context, the F flatness equations 
are, W = 0, Fi = dW/dcpi = 0, where W is built from holomorphic monomials invariant under the complete set of 
non-Abelian and Abelian gauge symmetries, including the global type string theory symmetries. 

Most of our applications deal with the fleld space of non-Abelian singlets charged under the Abelian gauge groups, 
G = {Y[a=i^O-)a) ^ U{1)a- The D-flat directions are parameterized by the monomials. Pa = Yii^'l^ ; solving the 
simultaneous equations: Qa{Pa) = 0, QAiPa) 7^ of sign opposite to ^a- By contrast, the superpotential couplings 
are constructed from the gauge invariant monomials with respect to both non-Abelian and Abelian (anomaly free 
and anomalous) gauge groups which obey the string selection rules, listed for the orbifold Zq-u in Eqs. (B.IO). 
The submanifold of vacua for non-Abelian singlets [55] is parameterized by the holomorphic invariant monomials. 
Pa = Yli 'Pi' ) b'i ^ 0] described by the column vectors, — (rf , • • • , r^) G Z!^, subject to the equations, 

BP 

Qa{Pa) = -QT:iQ-<t>)i = Pa{<t>) X ^rfQl = 0, ^^(Pa) = Pa{<t>) X J^^'i^A = 'CaU, [cA > 0]. (A.2) 

i i i 

These equations are solved by the simple solution, \(f>i\/^yrf — \(j>\ — 1/ca- The flat directions may have 0, 1 

or more dimensions, where the 0-d case involves a flxed common VEV, (j) = (— Ca/ X^j Oa''? the 1-d or 
higher dimension cases involve a single or more free complex parameters. When the system of linear equations is 

underdetermined, the solutions depend on free continuous parameters whose number identifles with the moduli space 
dimension, V < N — rank{A). The vector space of column vectors, r, is generated by the basis of column vectors, 

jfA _ ^^A^ . . . ^ ^A-j^ associated to the holomorphic invariant monomials. Ma = Di "^I' : = 1) " " ■ > 2-*] such that any 
solution can be written as, P" = Jla^^A^j for n > 1 and ua E Z (positive or negative signs). More conveniently, 
one can also consider the superbasis of one-dimensional invariant monomials. Pa, satisfying the condition that they 
cannot be factorized into products of two or more invariant monomials. 

The D flat direction described by the monomial solution. Pa {(j)) = Hi 4>i'^ > can be lifted by F-terms of type 

A, = (Iliep <^i")" type B, Wf = ^(Yl'i^p ^ Pa] consisting of gauge invariant superpotential 

monomials with all, or all but one, field factors included in (</>). In the terminology of [56, 57, 62], these refer to the 
gauge invariant non-renormalizable couplings in the field theory action, prior to applying the string selection rules. 
For the F flatness to remain valid to arbitrary orders of the full superpotential, an infinite number of conditions must 
be imposed for a type A direction but a finite number for a type B direction. The F fiatness conditions are, of course, 
more restrictive when set on the string theory action. 

Appendix B: Review of string construction for Zq-h orbifold 

The orbifold Zq-h is described by the rotation vector, 6^ = e'^^'^'"^ , or the twist vector, w| = i(l,2, — 3), acting 
on the complexified coordinate and fermion field components , tp^ associated to the complex planes of Tf. This 
orbifold is equivalent to the direct product of suborbifolds, Z2 x Z3, with twist vectors, w| = Swg = ^(1, 0, —1), w| = 
2vl = ^(1,-1,0). The anisotropic compactifications [10, 12] are realized with 6-d tori which factorize on the three 
maximal tori Tf of the rank 2 semi-simple groups, G2, SU{3), 50(4), with lattice basis vectors ei,2, 63,4, 65,6, and 
orbifold point group action represented by the 2x2 matrices for the Coxeter operators Cj c:i @^ in the lattice bases, 

0(G2) = (_\ ,^2) ' = (-1 -1) ' ©(^^W) = . (B.l) 



The input data for the orbifold fixed points and shift vectors in the various twisted sectors is displayed in Table VI. 
The shift vectors Ugj G A, associated to the fixed points / in sectors Tg ~ Q^, are defined by, (1 — 6^)/ = Ugj, modulo 
elements of the sublatticc, A^g = (1 — 6^)A, of the torus lattice A. The shift vectors arc in 1-to-l correspondence with 
the orbifold space group elements, {Q^,Ugj), which act on the orthogonal frame coordinates as, X — > Q^X + %,/. 
Each shift vector Ugj represents a conjugacy class element of the torus lattice coset, A/A^g. The fixed points /, shift 
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TABLE VI: Geometric data for the orbifold Ze-j/ on the torus = x Ti x T| of lattice, A(G2) + A{SU{3)) + A{SO{i)). 
The line entries refer to the twisted sectors, Tg, [5 = 1, ■ ■ ■ , 4] of Ze-ii with Ti ~ Ts, and the associated Tg-^^g^ sectors of the 
orbifold, x Za. The second colum displays the number of fixed points Mj^ in the three complex planes Tf . The third column 
consists of three subcolumns which display for the lattice A(G2), the fixed points and shift vectors, fa, Ua and fab, Uab, for 
the elements G^''' and Q^, the linear combination eigenvectors, 17(3) >, and the eigenvalues 7(5), with ui = e^"^'^^. The fourth 
column consists of two subcolumns which display for the lattice A(SU{3)), the fixed points and shift vectors, /„3, m„3, for 
the elements and G^'**, and the range of na. The fifth column consists of three subcolumns which display for the lattice 
A(S'0(4)), the fixed points and shift vectors, fn2,n'^, Wna.na' ^'^^ elements and G^, and the range of n2, n'2- The last 
column displays the modes multiplicities V assigned in the various twisted sectors Tg, in the case W2 = with unresolved 
quantum number 1x2=0, 1 for the orbifold group action on . The entries for T> correspond to the entries for the eigenvectors, 
17(g) >, of eigenvalues 7(3). 



Tl = A(G2) 



|7(S) > 



7(s) 



T| = A(S!7(3)) na 



= A(SO(4)) n2 n'2 



Tl,2+ri 1 



T'2+T4 
'^0,1 +'^0,2 



T3 
"^1,0 



1,3,4 



3,3, 1 



/a = f ei, ua = - 
[0 = 0,1,2] 



/„t = ^(aei+i,e2) 
Wjj(, = aei+6e2 
[o,!.=0,l] 



//0>;^(/l+/2)x / 1 N 



'/( 



/OI+/1O + /1I 

/oi+^'/oi+^' /ll 



_ [0. f e3 + 3<!4 

3=3 + §=4l 0,1,2 
^3=0,63,63+64 



0,1,2 



= -T-'5 + -T-'=e 



fc=0 
Ua=0 



0,1 0,1 







0,1 0,1 



vectors Uf and sublattices Agg = (1 — 6^)A(T|) of the twisted sectors g are defined as follows for the three 2-d tori, 
Tf, with G2, SU{S), SO{A) group weight lattices: 

A(G2) : (1 - e^'^)fa ^ua^ -ae2, fa = ^ei, [a = 0, 1, 2]; (1 - e^)/^,^ = w^,^ = |ei + ^62, [a, h = (0, 1)]; 
Ae ^ (ei, 62), A^2,4 ~ (ei, 862), Aqs ~ (2ei, 2e2). 

A(^C/(3)): (l-e^)/„3=«„3 = [e3,e3 + e4]=n3e3, /„3 = [0,^^^^, ^^^^], [713 = 0,1,2]; 
A|g - (63-64,863), [g= 1,2,4]. 

A(5'0(4)) : (1 - &^)f„2,„'2 = "«2,ni = -"265 - 71266, fn2,n'2 = 2 ('^265 + n'aee), = 0, 1]; 

A|, -(265,266), b= 1,3]. (B.2) 

The orbifold fundamental group allows for three discrete Wilson lines: W2, W2 of order N2 = N2 = 2, around two 
dual onc-cyclcs of Tg , and of order = 3 around the 63 one-cycle of torus T2. Upon turning on the Wilson lines, 
the Tg twisted sectors split into twisted subsectors, (3,7,7x3,712,712), labelled by the discrete parameters, 713 =0,1,2 
and 7i2 = 0, 1, n'a = 0, 1, for the lattices A(T2) ~ A(5J7(3)) and A(T3) ~ A(SO(4)). The complex phase parameters, 
7(3), describing the action of on the fixed points of the lattice A(T^) ^ A(G2) are needed to specify the degeneracy 
of string modes in the sectors <? = 2, 4 and <? = 3. b The string states arc described by the 50(8) group weight vectors, 
r°, [a = 1, • • • , 4] for the right-moving fermion fields; the Es x group weight vectors, , [I = 1, • • • , 16] for the 
left-moving 16 coordinate fields; and the oscillator numbers, N.^-, Nf" & Z, [i = 1,1; 1,1= 1,2,8] for the 8 right- 
and left-moving complex coordinate fields of T^. The string squared mass spectrum is evaluated in terms of these 
quantum numbers by means of the formula 

^ = ^(iv,M^) + 4^)) + i ^(r'^ + gv^r - 4'^ - 1, 

* / a=l 

[-^0^^ " ^ X] - \9Vi\), gvi = gvi mod l = gvi- [gvi], 
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TABLE VII: The SO{8) group weight vectors, ri''^ and rl'^, assigned to the massless right moving modes, M% = 0, with left 

and right chirahtics (/, r). The vector and spinor representations, (v, s) differ by the weight vectors pi^r assigned to the 
supercharge generators, ri''' = r'y^ + pi,r, \pi,r = (±| ± | ± |, ±|)]- The line entries refer to the untwisted and twisted sectors, 
g = 0,1, ■ ■ ■ ,4. No massless mode solutions arise for right movers carrying oscillator exitations, Np' j. 



g r[ r' = r„ + pt 


g = r,, + pjj 


(100,0) (+1-1-1,-1) 


(-100,0) (-i + i + i,+i) 


1 (001,0) (-1-1 + 1,-1) 

2 (001,0) (-i-l + i,-i) 

3 (0-12,0) (-1-I + I,-!) 

4 (0-12,0) (-1-1 + 1,-1) 


5 (-1-22,0) (-1-1 + 1, +1) 

2 (-1-11,0) (-l-l|,+i) 

3 (-1-11,0) (-i-i|,+l) 

4 (-1-22,0) (-1-2 + f,+l) 



Xg,n, =9V + Hf^^Wa = gV + nsWs + n^W^ + n'^W^] (B.3) 
(in units of mg) where Xg^nf denote the shift vectors in the Eg x Eg group lattice depending on the orbifold and Wilson 
lines gauge cmbeddings, V, W3, W2, The oscillator energies are defined by ujf^ = gvj for gvj > and 1 — [(/u/j for 
gvj < 0, with N^'^ G Z+. An equivalent definition is, uj^f^ = gvj mod 1, Wj^'* = —gvi mod 1, for the determinations 

obeying, < Lo[^f < 1. In the Zq-u orbifold, the string states generally occur in CPT conjugate pairs of opposite 4-d 

chirality for the sectors, g and g' = 6 — g, except for the twisted sectors T2 and T4 with g = 2, 4, which admit both 
left and right chirality modes each (because of the supcrsymmctry A/" = 2), and the (self-conjugate) twisted sector T3 
with 3 = 3. The vector and spinor weight vectors r^g of the SO{8) symmetry group for right moving fermions are 
associated to boson and fermion superpartners of the massless chiral supermultiplets as displayed in Table VII. 



The GSO projection on the orbifold singlet states is determined via the one-loop partition function by the condition, 
P{g,nf,^,(t)) = 1, where 

1 

P{9,nf,l,(t>) = ]y D A''(3,n/,7,(/)), 

[/^^{g,nf,^,cl))=^{g,h)<i){g,h)e^'''^^''+''<>'-t^-'''''"s-^''+^^^^^^ (B.4) 

The twists along the world sheet spatial and temporal directions are denoted by g and h and the terms A'' include the 
complex phases, ^{g,h) = 7'' (5), eigenvalues of the orbifold twist action on the fixed points in the 2-d torus lattice 
A(Tf ) = A((S'2), and the complex phases from oscillator excitations, (f){g, h) = 4''^{g), 

(j){g) = e''''^^='.J^^'~^^^^\ [4>i =visgn{gvi), $j ^ -visgn{gvi), /, / = 1, • • • , 3]. (B.5) 

The level matching for right and left movers entails the conditions, N{Xg^^^ — (gv)"^)) = mod 2, with additional 
conditions involving the scalar products of Wilson lines, W3, W2, W^- For the twisted sectors, T2,4 and T3, the 
projections for the and Z2 suborbifolds can be implemented by summing over the time twists in the subsectors 
{h = 3, TO2, m'2) and {h = 2, ms) as follows 



3=(2, 4), /l = 3: P{g,n^,^,<i))=^^{g)<l/'{g)^ ^ ^2^^[{P+U9V+n,W,))■(hV+n^^W^+rn'^W^)-(r+^av)■hv\ ^ 

1712 rfn^ —0.1 

g = 3,h = 2: P{g = 3,n2,n'2,^, 4>) = j''{g)ct>\g)\ V e^'''^^''+i^BV+n,w,+n',wi)).{hv+msWs)-{r+yv).hv] 



m3=0,l,2 

(B.6) 

Except for the singly twisted sector, Ti, where the conditions in Eq. (B.4) hold with fixed values of ns, n2, the 
projections in the other sectors can be concisely stated in terms of simple sets of conditions on the weight vectors 
P^. For the untwisted sector, these conditions are: P -V € Z, P ■W2€ Z, P • € P • W3 e For the twisted 
sectors, T2^4 and T^, the projections are more easily implemented by imposing the conditions 



(2, 4), h = 3: 



hP-W2GZ, hP-W^GZ, [P = P + Xg,„3,o,o] 
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• 5 = 3, /i = 2: ^'^^ . , , . (B.7) 

V /iP-VFse^, [P = P + Xg,o,„2,n^] / 

This is the prescription for the GSO orbifold phase that we have used to obtain the string spectra of models A\, B. 

Use of the world sheet modular invariancc (level matching conditions) allows replacing the conditions on the scalar 
products in Eqs. (B.7) by the simpler ones [10]: P • VF2 G Z, P • e Z for 5 = (2, 4) and P • W3 G Z for g = 3. The 
following alternative projections are also proposed in [15] 



• 5 = (2, 4) 



us are aiso proposeu in [loj 

/3[(7V/' - Nf)vi + q^{g) - {r + gv) ■ v + {P + gV + nsWs) ■V]€Z\ 
[ ^p + gV + nsWs)-{Z'^€Z J' 

= 3- (Vi'-Ni')vi + q^i9)-{r + gv)-v + {P + gV + n,W2+n'^W^)-V]eZ\ 
*^ ■ V {P + gV + n2W2 + n'^W!^)-Wz&Z J' ^ 

where j{g) = e^"'''<^^\ It is important to remark at this point that the orbifold projection for sectors with fixed 
subtori boils down to the invariance tuider the world sheet modular group. The modular invariance is vital if the 
mixed gauge and gravitational anomalies in the various Abelian gauge factors U{l)a are to satisiy the universality 
relations [81, 84] 

\j^Tr{Q^) = -^Tr{Ql) = ^Tr{QMR)) = -^Tr{QaQl) = 8nH^GS, (B.9) 



24/e;(j 3/ua 2/ua ka kit 

which are necessary for the Green-Schwarz type anomaly cancellation, mechanism to work. When the orbifold space 

group action by h = {Q^, Uhj) along the string time direction docs not commute with that of the constructing orbifold 
element along the string space direction, g — (O^ ,Ugj), the physical states must be constructed by summing over 
elements of the orbit, h^gh~'^, and including suitable complex phase factors to ensure the SL{2, Z) modular group 
invariance [17, 20, 85]. This situation indeed occurs for the sectors, 5 = 2,4 and 5 = 3, where the fixed points in Tl 
with lattice A(G2) are reshuffled by the action of h. The physical string modes decompose then into eigenvectors of 
es, [g = 2,3,4] with eigenvalues 7(5) = ±1 in the [g = 2, 4) sectors T2, T4 and 7(5) = I.oj.uj"^, [uo = e^^'^/^] in the 
(5 = 3) sector T3. The degeneracy of eigenstates is 2? = 2 in Ti; 2? = (1, 2) for 7 = (1, -1) in T2, T4,, and V = (4, 2, 2) 
for 7 = (l,w^,u;) in Ta, as displayed in Table VI. 

For non-commuting twists, [h,g\ 7^ 0, the complex phase "f{g,h) is not determined in a unique way by the GSO 
projection. Since 7(5, h) fixes the degeneracies of string modes, its choice has clearly an incidence on the U{l)a 
groups anomalies and hence on the modular invariance. The choice of 7(5, h) can be set uniquely by requiring the 
mixed gauge and gravitational anomalies to satisfy the anomalies universality relations in Eq. (B.9). Given a massless 
string state, specified by the quantum numbers, P^, , TV/", then the freedom on j{g,h) is taken into account by 
redefining, 7(9, h) x{^)l{9: in terms of the GSO modified condition, A''x'' = 1, so that the mode multiplicities 
determined by 7(5, h) give anomaly coefficients satisfying the requisite universality relations. Although the natural 
choice, X = 1, applies for the large majority of modes, there do occur cases where a non- trivial x in is needed to 
satisfy the universality relations. The factors x{h) depend on the prescription used for the modified gamma phases. 
With our GSO projection prescrition in Eqs. (B.7), one exception occurs for three modes of sector in model A2 and 
another one for two modes of sectors T2 and T4 in model B. 

Finally, we quote the string selection rules on the n-point world sheet correlators of vertex operators, < 
niLi ^(-^/i ^i) >• Applied to the superpotential couplings, the conditions set by the gauge symmetries, the H- 
momentum conservation, and the orbifold point and space symmetry groups read: 



^[P^(0 + X^',„,(0] e AiEs X Es), [Xgj =gV + n^jWa]. 
1=1 

n 

^Ri{l) - 1] = mod ( — ) = mod (6,3,2), [R\1) = r\l) + giv\l) - N^{1) + N^{1)\. 
(=1 l^^l 

n n 

= mod 1 =^ "Y^gi = mod N. 

1=1 1=1 

n n n 

^(l-ef')(/z-A/)=0modiV =^ ^713(0 = mod 3, ^(n2(0, "2(0) = (0> 0) mod 2. (B.IO) 



1=1 



Some of the above selection rules can be formulated as ordinary [11, 86] or R type [87] discrete symmetries. We have 

omitted from the list in Eq.(B.lO) the constraints on couplings from the space orbifold group action on the torus 
with lattice A(G2). These consist of the following two selection rules, stated around Eq. (C.9) of [10]: 
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• Rule / on the gamma phases of physical modes: HlLi ll = 1- Using the definition, 7; = 7(5;) = e^™'^'<i , this can 
also be written as, X^ILi = mod 1. 

• Rule // for couplings involving the twisted sector modes T2,4 only, or twisted sector modes T3 only, times 
untwisted sector states. The version of this rule derived in [13] states that the column vector of gamma 
phases for fields in n-point couplings T2\T^^ or T^Tq^, with pi 7^ 0, n > p2 > must obey the condition, 
Qii ^ [(P) 0) • • • ) 0) + perms], [p ^ 0]. Thus, the allowed couplings are those containing at least two physical 
modes with non-trivial gamma phases, ji ^ 1. The version of this rule derived in [10] uses instead the conjugacy 
classes of the lattice cosets, A/Ae?. associated to the fixed points /^^ = yei, [a; = 0, 1,2] in sectors g = 2,4 
and /oj^fcj = ^ei + ^62, [a; = 0, 1; bi = 0, 1] in sector g = 3, in the notational conventions of Table VI and 
Eq. (B.2). In terms of the fixed point content of the physical states (with fixed 7; eigenvalues), the orbifold space 
group constraints arc described for the couplings ^ by the Z3 group selection rule, 5^"=! '^■i = *-* mod 3, 
and for the couplings T^^Tq^ by the Z2 x Z2 group selection rule, X^JLi = mod 2, Yl^=i = niod 2. 

Convincing arguments are given in [13, 20, 85] that Rule / is not a genuine selection rule, since it is found to be 
automatically satisfied by physical states obeying the GSO orbifold projection, ri/, 7, 0) ~ 1, once the selection 
rules on P^{1). (l) are imposed. In a consistent string model complying with the orbifold GSO projection, 

and hence with the world sheet modular group invariance. Rule / is redundant. However, the redundancy was proved 
with the specific prescription for the orbifold projection quoted in Eq.(B.8). Since our prescription in Eq.(B.7) is 
quite close but not strictly identical to that used in [13, 20], there is no guarantee that Rule / is automatically implied 
by other rules in our case too. That a selection rule is superfluous in one projection prescription and not in another 
appears odd at first sight but cannot be excluded. It certainly is a logical possibility if several consistent string models 
arise from the same data set of shift vectors, V, W2, W3. To settle this issue on practical grormds, we have applied 
our search procedure for models Al, B with the rule OILi 7' — 1 included, and fomid that this did forbid certain 
superpotential couplings allowed by the rules in Eq.(B.lO). However, the excluded couplings represent a very tiny 
fraction of the large set of allowed couplings, which cause insignificant changes on the effective couplings. We have 
also tested the Rule II above and found that this had a negligible incidence on the allowed couplings. 
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